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Abstract 



Recently, there has been growing interest on primordial "blue" (n > 1) perturbation spectra, 
motivated both by a composite set of observational data on large scales and from the point of 
view of theoretical model building. After reviewing the theoretical (inflationary) motivations 
for these blue spectra, we consider various observational constraints, within both Cold Dark 
Matter and Mixed Dark Matter scenarios. In particular, using linear theory, we discuss large- 
scale bulk flows and the X-ray cluster abundance. We also perform various N-body simulations 
of these models to study the clustering properties of the matter distribution and the peculiar 
velocity held. 

1 Introduction 

The anisotropy of the Cosmic Microwave Background (CMB), detected on large angular scales 
by COBE (Smoot et al. 1992; Bennett et al. 1994), provides an insight on primordial pertur- 
bations. This measurement yields, on one hand, the amplitude of density perturbations at the 
largest observable scales, and, on the other hand, the spectral index n, that fully determines the 
power-spectrum [i.e. P{k) oc Ak^] of the primordial density fluctuation fleld, up to a possible 
contribution by a stochastic background of gravitational waves. These constraints, promptly 
taken into account in the current cosmological debate, can discriminate among the possible 
scenarios. In particular, the standard (i.e. n = 1) Cold Dark Matter (hereafter SCDM) model 
with the COBE normalization predicts too much structure on scales ^ 10 Mpc {h being 
the Hubble constant in units of 100 km s^^ Mpc^^), e.g. in terms of an excessive abundance 
of clusters (e.g. White, Efstathiou, & Frenk 1993; Colafrancesco & Vittorio 1994), and of too 
a large rms pairwise galaxy velocities on Mpc scales (Davis et al. 1985; see however Mo, Jing, 
& Borner 1993; Zurek et al. 1994). Tilted (i.e. n < 1) CDM models are in better agreement 
with observations on small scales, but the amount of tilting is strongly constrained from data 
on the large-scale velocity fleld (Moscardini et al. 1995) and by cluster clustering properties 
(Phonis et al. 1995). The most severe flaw of tilted CDM models is the very late epoch of 
galaxy formation. 

A fashionable alternative to SCDM, recently emerged to overcome these problems, is pro- 
vided by Mixed (i.e. cold + hot) Dark Matter (MDM) models: adding to the cold dark matter 
a small contribution (~ 10 — 30% of the critical density) of massive neutrinos erases, via free- 
streaming, part of the power on the scales of galaxies and clusters. This scenario will also 
be reinforced from the particle physics side, if recent claims for a mass of order 2 eV for the 
ji neutrino will be conflrmed (Louis 1994). The standard MDM model is based on the same 
assumptions as SCDM: in particular, the primordial power-spectrum is assumed again to be 
scale-invariant (i.e. n = 1). Although the results of linear analyses and numerical simulations 
(e.g. Klypin et al. 1993; see also Davis, Summers, & Schlegel 1992; Jing et al. 1994; Cen & 
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Ostriker 1994; Klypin, Nolthenius, & Primack 1995) are certainly promising, some problems are 
clearly present: an example is the excessive delay of the epoch of galaxy formation, in contrast 
with observations of quasars and damped Lyct system at redshifts z > 3 (Ma & Bertschinger 
1994; Klypin et al. 1995). 

The analysis of the first two years of the COBE DMR data indicates a most likely value 
for the spectral index in the range n = 0.9 — 1.9 (Bennett ct al. 1994; see also Gorski et al. 
1994), a value also suggested by the Tenerife experiment at 5° (Hancock et al. 1994) and by 
the MAX measurements at 1° (e.g. Gundersen et al. 1993; Devlin et al. 1994). An n > 1 
perturbation spectrum can also be at the origin of the claimed bulk flow on scales of 150 
Mpc (Lauer & Postman 1994) and of the large voids in the galaxy distribution of the CfA 
survey (de Lapparent, Gcllcr, & Huchra 1986; Geller & Huchra 1989), at scales of order 50 
Mpc, as pointed out by Piran ct al. (1993). On the other hand, it is also possible to put 
upper bounds on the value of n. Carr & Lidsey (1993) found n < 1.5, from primordial black 
hole overproduction constraints on nucleosynthesis and the COBE data (see also Garr, Gilbert 
& Lidsey 1994). Hu, Scott, & Silk (1994) set an independent constraint, n < 1.54, by using 
recent experimental limits on the CMB spectral distortions found by COBE FIRAS. 

Most theories for the origin of the primordial fluctuations predict a spectral index n <1. In 
particular, inflationary models based on simple inflaton potentials originate this type of spectra. 
Nevertheless, in recent years some authors have pointed out that suitable inflationary dynamics 
can easily account for blue spectra (Liddle & Lyth 1993; Linde 1994; Mollerach, Matarrese, & 
Lucchin 1994; Gopeland et al. 1994), in particular, in the framework of the so-called hybrid 
models. 

In this paper we investigate the predictions of a possible variant of MDM models, by assum- 
ing a "blue", i.e. n > 1, primordial density perturbation spectrum; this BMDM model (where 
"B" now stands for "blue"), suitably normahzed to COBE DMR, may alleviate the small-scale 
problem of the standard (i.e. n — 1) MDM model. An important advantage of "anti-tilting" 
the spectrum is that it provides an earlier galaxy formation epoch. Here we restrict ourselves to 
discuss the present matter and light distribution, both in the linear and nonlinear regime. We 
refer to de Bernardis, Baldi, & Vittorio (1995) for a discussion of how MDM models compare 
with the anisotropy data at degree angular scales. 

The plan of the paper is as follows. In Section 2 we discuss the theoretical motivation for 
blue power-spectra from the inflationary point of view. In Section 3 we show the predictions 
of the linear theory for the large-scale peculiar velocities and for the galaxy cluster abundance. 
In Section 4 we present the analysis of our N-body simulations. Finally, in Section 5 we discuss 
our results and draw our main conclusions. 
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2 Inflationary models 



Quantum fluctuations of scalar flelds during inflation provide a very appealing explanation for 
the origin of the primordial density perturbations. It is useful to describe these perturbations 
in terms of the gauge-invariant variable Its power-spectrum to first order in the slow roll 
approximation is given by 



(1) 



aH=k 



1 /2 

Its spectral index, defined by a = dlnP^' /dink, is related to the spectral index of density 
perturbations on constant time hypersurfaces before matter-radiation equality as n = 2q; + 1. 
Most infiationary models predict a spectral index n < 1. This is the case for example for models 
based on a polynomial potential, which give rise to a power-spectrum P^^'^{k) logarithmically 
decreasing with k. Another example is the power-law infiation, which is based on an exponential 
potential and gives rise to a power-spectrum decreasing with k as a power (Lucchin & Matarrese 
1985). This is due to the fact that in most infiationary models decreases with time as 

the infiaton evolves, thus the larger wavenumber perturbations, that left the horizon later, have 
smaller amplitude. However, as it has recently been discussed by MoUerach, Matarrese, & 
Lucchin (1994), there is a class of simple inflationary models leading to density perturbations 
with n > 1 (see also Liddle & Lyth 1993; Carr & Lidsey 1993; Linde 1994; Copeland et al. 
1994). These are produced while the inflaton field rolls down to a potential minimum with 
potential energy dominated by a cosmological constant-like piece. In this case \<p\ decreases as 
(f) approaches the potential minimum, while remains approximately constant in such a way 
that H^/\4>\ increases with time. The shape of the infiaton potential leading to P^^'^{k) oc k"' 
with a = const > has been obtained by MoUerach, Matarrese, & Lucchin (1994); it looks 
quite complicated. However, if the condition that a is constant for all scales is relaxed, a huge 
class of potentials with a non-vanishing energy density in the minimum, which corresponds to 
an effective cosmological constant, can work. For example, the potentials V{(f)) — Vq + 
give rise to a spectrum decreasing with k for wavelengths that left the Hubble radius when 

where the right hand side has to be evaluated at the time when the wavelengths of interest left 
the horizon during infiation. In particular, we see that the quadratic potential case, m — 2, 
gives rise to an approximately constant spectral index n > 1. 

One can think that these kinds of potentials will have problems to end infiation and reheat 
the universe, as the two usual reheating mechanisms, namely rapid oscillations of the infiaton 
field around the potential minimum and first-order phase transitions with bubble production 
do not work in this case. This reheating problem can nevertheless be avoided by the "hybrid 
infiation" mechanism, recently proposed by Linde (1991a,b, 1994). According to it, infiation 
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ends in the following way: the slow roll of the inflaton field at a given time triggers a second- 
order phase transition of a second field, whose false vacuum energy density is responsible for 
the cosmological constant term. The field rolls down to its potential minimum and oscillates 
around it. A particular potential in which this scenario can be realized is (Linde 1991a,b, 1994) 

a) = ^(M^ - + + (3) 

For (f) > (pc ^ M/g, a is in its false vacuum at cr = 0. The potential for in this regime is 
given by V^(0) = M^/4A + m^0^/2. As we have discussed, when the cosmological constant term 
dominates, this gives rise to an inflationary phase with a power-spectrum Pc^{k) increasing with 
k. From eq.(0), the spectral index for the perturbations is given by a ~ AXni^ / . To be 
consistent with the approximations done, this result holds only for small values of a. For a 
given value of a, this fixes the relation between and M^. When (p < (pc, cr rolls down to its 
minimum and the universe reheats. 

There are some constraints on the model parameters in order that it works. In order to 
have at least 60 e-foldings of inflation for > 0c, one needs ^ rn^Xe^^^Mp/ g, and in order 
that inflation finishes soon after reaches 0c, one needs ^ 3A^Mpm^/27r^. 

The amplitude of the perturbations that left the horizon around 60 e-foldings before the 
end of inflation (the ones relevant for structure formation) is given by 

We normalize its amplitude using the COBE Qrms-ps determination, P^(/c6o) = A^{n) = 
(Qrms-Ps/To)Mor2(2 - n/2)r(9/2 - n/2)/7rr(3 - n)r(3/2 + n/2). From the analysis of the 
two first years of COBE data, Qrms-ps = 18.2e°-^s(^-"VK (Bennett et al. 1994) and the FIRAS 
value To = 2.726 ±0.010 K (Mather et al. 1994), the resulting values of A{n) are of the order of 
4 X 10~^ for the values of n of interest. This normalization, combined with the previous relation 
between and in terms of a, gives the values of m and M: 

M ~ 4.3v4(n)ae6°"— Mp, m - A7A\n)a^/^e'^'^''^Mp. (5) 
9 9 

For a given value of a, these relations fix the values of m and M. We have to check that the 
constraints imposed on the model parameters are satisfied. The condition M^0^/2 < M^/4X 
translates to A^/^ < 50e-^°°g/a^/^. We can ensure that this is satisfied for all the values of a 
of interest taking for the coupling constants g ~ C'(l) and A small enough. The condition for 
inflation to end soon after ~ 0c holds only for a ^ 0.17 (n ^ 1.35). The additional constraint 
that the energy density be smaller that the Planck density, M^/4A < Mp, always holds when 
the above conditions are satisfied. 
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3 Linear Tests of Blue Primordial Spectra 



Following the inflationary paradigm, we restrict ourselves to flat cosmological models, where a 
Hubble parameter h = 0.5 is assumed and baryons contribute only for the 5% to the critical 
density (i.e. fi^ = 0.05). The fraction of the critical density in cold and hot dark matter is 
^CDM = l — flB — ^u and fl^, respectively. The family of flat MDM model is then characterized 
only by the parameter Q^, that in this Section we let vary between (pure CDM) and ^ 0.5 
(the MDM model where hierarchical clustering still occurs). 

In order to study the present matter distribution in MDM models, we evaluate the transfer 
function of density fluctuations in each of the four components (CDM, HDM, baryons, photons), 
by numerically integrating the full set of equations for the perturbed matter-energy density fleld 
in the synchronous gauge formalism. The details of our calculation are described in de Gasperis, 
Muciaccia, & Vittorio (1995; see also Peebles 1981; Bonometto et al. 1984). The results for 
Jlj, = 0.3 are shown in Figure 1, where we plot the transfer function T{k) for each separate 
component: CDM, HDM and baryons, respectively. Here we assume that there is only one 
family of massive neutrinos. So, — 0.3 corresponds to a neutrino mass rrii, ^ 7eV. We also 
plot the total transfer function defined as follows: 

TTorik) = flcDuTcDMik) + Vt^T^{k) + VtBTB{k) . (6) 

The transfer function is given in arbitrary units. However the relative amplitude of T{k) 
at different redshifts reflects the linear growth of density fluctuation in an Einstein-de Sitter 
universe. We evaluated the transfer functions at redshift z — 20, in order to properly start our 
N-body calculations (see next section), with the right statistical segregation among the different 
components. In Figure 1 it is evident that baryons (dotted line) catch up completely the CDM 
(dashed line) at the present time. The transfer function of massive neutrinos is initially lower 
than the CDM one because of free-streaming. However, it gets presently comparable to that 
of the other components, at least for the range of wavenumbers shown in Figure 1. 

Figure 2 shows the rms value of the density fluctuation fleld at the present time, a{R), 
smoothed with a top-hat fllter of radius R. We plot (t{E) for the following models: SCDM, 
a blue (with n = 1.2) CDM model (hereafter BCDM), MDM and BMDM (n = 1.2) with 

= 0.3. All models are normalized to the rms value (30 fiK) of the temperature fluctuations 
in the DMR maps. It is apparent the reduction of power on small scales and the flattening 
of o'(-R) for the MDM model compared to SCDM. Our BMDM model is intermediate between 
SCDM and MDM. So, this model can in principle be consistent with the small-scale clustering of 
galaxies. This alleviates the problem of the standard MDM model, where structure formation is 
expected to occur later than in SCDM. To roughly quantify this issue we evaluate the redshift of 
formation of structures by asking when, on a given scale, a one sigma fluctuation goes nonlinear. 
In Table 1 we show the redshifts of nonlinearity for three mass scales (M = 10^^, 10^^, and 
10^^ M0), corresponding to typical masses of galaxies {zni,gai), groups {zni,gr) and clusters of 
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galaxies {zni,d), respectively. 

Table 1. Redshift of Galaxy, Group and Cluster Formation 





^nl,gal 


^nl,gr 


^nl,cl 


SCDM 


4.4 


2.5 


1.0 


BCDM 


11.4 


6.5 


3.1 


MDM 


0.6 


0.3 


0.1 


BMDM 


1.9 


1.3 


0.7 



3.1 Bulk Velocities 

We test our models against peculiar velocity data. We consider the Local Group (LG) velocity 
relative to the CMB (Kogut et al. 1993), assigning a velocity of \v^^^\ = 627 ± 22 km s"^ to 
a sphere of radius Ri = 5 Mpc. We also use recent data by Dckcl (1994), which give the 
bulk flow (m = 2, 7) of spheres of radius Rm = 10, 20, . . . , 60 Mpc. These data 

were derived after reconstructing with the POTENT method (Bertschinger et al. 1990) the 
three-dimensional velocity field smoothed by a Gaussian window with a radius of 12 Mpc. 
Therefore, in the following we will consider the set {v^"^^} (m = 1, 7) corresponding to the LG 
velocity and six bulk flows. 

The likehhood of the data p({'y^"^)}|7i) under the hypothesis H, in our case a model with a 
given ill, and n, can be written as 



(27r)3^/2 1 1/ 1 3/2 



exp 



AT N 



^ i=i j=i 



(7) 



where Vj~^ and \V\ are the inverse and the determinant of the 7x7 correlation matrix, Vij = 
^^») . vO)y^^ respectively. The condition H is contained in the elements of the correlation 
matrix, which can be written as 



= dkP{k)Wi{kRi)Wj{kRj) + ^Sij 



(8) 



where the indices i and j refer to the velocities of spheres of different radii as indicated in Table 
2 and P{k) is the density fluctuation power-spectrum of the chosen theoretical model (properly 
normalized to CORE DMR data). Moreover, in eq.(8), W,{kRi) = 3ji{kRi)/{kRif. 

The observational data \v^^s\ with their error bars Hg^g are compared in Table 2 with the 

theoretical predictions for the three-dimensional rms bulk flows \vj-}}\ [given by \/3Vii {i = 1, 7)] 
obtained for SCDM, BCDM (with n = 1.2), MDM and BMDM (n = 1.2) with = 0.3. The 
corresponding values of the bias parameter b = (T~^(8/i~^Mpc) are also given in Table 2. 
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Table 2. Bulk velocities. Observational data and theoretical predictions. 











SCDM 


BCDM 


MDM 


BMDM 








bias 


0.95 


0.50 


1.47 


0.91 


i 


Ri 


1 1 


y{i) 
^obs 


i-Wi 


-<i) 














1 


5 


627 


22 


880 


1387 


723 


1107 


2 


10 


494 


170 


509 


705 


485 




670 


3 


20 


475 


160 


467 


63J 




452 


614 


4 


30 


413 


150 


418 


562 


409 


545 




5 


40 


369 


150 


377 


496 


371 


48S 




6 


50 


325 


140 


338 


438 


335 


434 


7 


60 


300 


140 


310 


394 


308 


391 



Now we consider models with 1 < n < 1.5 and Q„ < 0.5 and we look for the model most 
consistent with the data, i.e. the one that maximizes the Likehhood of the set {v^^^^ohs- 
find that this happens for ft,^ — and n = 1.1, if we maximize the joint probability of having the 
measured bulk flows ano? the CMB dipole. We find instead fi^ = and n — 1.2, if we maximize 
the conditional probability of having the measured bulk flows, under the condition that the 
CMB dipole has its measured value. The conditional probability is obtained by dividing eq.(7) 

by 



1 Iv^^M^i 



(9) 



2 Vn 

Thus, BCDM models seem to be favoured from this kind of analysis. 

In order to define a 95% confidence level for and n we proceed as follows. We define a 
likelihood ratio A = p{{v^'^^}\'Hi)/p{{v^"^^}\'H2): and we test the condition Hi (a model with 
given 0,1, and n) against the condition 7i2 (the most probable model). Here p is either the joint 
[of Eq.(7)] or the conditional probabihty of reproducing the observed set {v^^^}f,i,s- We generate 
10,000 realizations of the set {v^"^^ under Hi and 'H2, respectively, and from them we evaluate 
p(A|7ii) and p(A|7Y2) (see Del Grande & Vittorio 1992 for more details on this method). We 
find that BMDM models, with a small fraction of hot dark matter, Q,y ^ 0.3, can provide a 
reasonable agreement with the peculiar velocity data. We will compare this conclusion with 
the analysis of detailed N-body simulations in Section 4.2. 



3.2 Cluster Abundance: the X— Ray Constraints 

In this section we test the predictions of models with primordial blue perturbation spectra 
against the local X-Ray Luminosity Function (XRLF) of galaxy clusters (Kowalski et al. 1984). 
In hierarchical scenarios the XRLF can be derived from the cluster mass distribution using the 
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appropriate M/L relation between the cluster mass M and its X-ray luminosity L. Several 
evidences (see e.g. Colafrancesco & Vittorio 1994 for a discussion) suggest that the cluster mass 
distribution can be reasonably represented by the Press & Schechter (1974) formula 



dM 



(27r)V2 



e-^'/2 ^ ^10) 



where Vc = Sc/<^{M, z), Sc = 1.7, and p is the mean density at time t (the mass M is defined 
through the corresponding radius of a top-hat filter of radius R) . The normalization is uncertain 
within a factor 1 ^ X ^ 2 (see e.g. Bond et al. 1991). A direct comparison of the mass 
distribution in eq.(lO) with the data of Bahcall & Cen (1992) is not straightforward, as the 
mass M in the Press & Schechter theory can be related to the estimated mass within a fixed 
radius, M4 = M(r < 1.5 Mpc), only if detailed information about the cluster density 
profiles is available (e.g. Antonuccio-Delogu & Colafrancesco 1994, for a discussion). On the 
contrary, the X-ray luminosities are mostly sensitive to the inner parts of the density profiles, 
as the specific emissivities [mainly due to thermal bremmstrahlung from the hot intergalactic 
medium (IGM)] are Cj, ~ y^°' exp{—y) , where y = hu/kT and a 0.3 — 0.4 (Gioia et al. 1990). 
So, apart from the theoretical difficulties in explaining the existence of a cluster core, the X-ray 
luminosity oc J drr'^ff{r)T^^'^{r) is a more rehable quantity than Ma, due to its relation to 
the total cluster mass M. Here n is the density of the IGM at temperature T and the integral 
is extended over the cluster volume. We adopt here the following parameterization for the local 
cluster X-ray luminosity 

L (X M^p'^IC{AE,z) , (11) 

where c = 4/3, d = 7/6, and IC{AE,z) = Jy^^^ dy ey{y)/ dy ej,{y) is the K-correction in 
the energy band Ay. This correction, applied to the HEAOl energy band 2 — 6 keV, yields, at 
2; ^ 0, a reduction of the total luminosities by a factor ~ 3 (see Colafrancesco & Vittorio 1994, 
for a more complete discussion). 

Given N{M, z) and the ratio M/L oc M^^'^ from eq.(ll), we obtain the comoving XRLF 

Ar(L, z) ^ N(M, z)^ = l^N[M(L, z, AE), z] . (12) 

We plot in Figure 3 the local XRLF predicted by our models, normalized to COBE DMR, 
choosing the canonical value 5c ~ 1.7, to select collapsed objects. The XRLFs predicted in all 
our models are well above the observational data, indicating that COBE normalized models 
are unable to reproduce the cluster abundance, quite independently of the value of the spectral 
index n. If we allow for 5c values larger than 1.7, we can reproduce the observed XRLF at 
a better level (see Colafrancesco & Vittorio 1994). In fact, the mass distribution in eq.(12) 
depends upon the combination 5ch and X. So for a given model (i.e. for a given n^, and n) we 
may derive 5ch and X by a direct fit to the XRLF data. The results of a Chi-square analysis 
are shown in Table 3. The theoretical XRLFs, obtained with the values of Table 3 for 5ch and 
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X are also shown in Figure 3. Note that for — 0.3 and n — \ the difference between the two 
approaches is quite small. 

A discussion of the cluster abundance obtained in N-body simulations will be presented in 
Section 4.3. 

Table 3. Parameters of the XRLF fit 



iS^v, n) 


Scb 


J 


A,mm 


(0.1,1.3) 


2.7 


1.3 


1.93 


(0.2,1.4) 


2.7 


1.3 


1.74 


(0.3,1.4) 


2.8 


1.7 


1.59 



4 N— body simulations 

To study the large-scale structure of the Universe in models with blue primordial perturbation 
spectra we ran N-body simulations of the matter distribution. We used a particle-mesh code 
with a box of side 260 Mpc (we adopt here h — 0.5). We considered three different models: 

standard CDM (SCDM), a blue CDM model (BCDM) with primordial spectral index n = 1.2 
and a blue Mixed Dark Matter model (BMDM), still with primordial spectral index n = 1.2, but 
with QcDM = 0.65, Qb = 0.05 and fij, = 0.3. We choose this last model, even if only marginally 
consistent with the linear analysis of the velocity field presented in Section 3.1, in order to study 
an extreme case where the effects of the hot component can be largely appreciated; moreover 
this choice allows a direct comparison with the results of Klypin et al. (1993) who simulated a 
model with the same but with a primordial spectral index n — 1. 

We define the present time by fixing the normalization to the value implied by COBE 
DMR data in the absence of any relevant gravitational wave contribution: this corresponds to 
approximately 6 = 1, 6 = 0.5 and h = 0.9, for SCDM, BCDM and BMDM respectively 

In the case of CDM models (SCDM and BCDM), we used 128^ grid-points and 128''^ 
particles. Initial conditions were set using the standard approach based on applying the 
Zel'dovich approximation and assuming Gaussian primordial fluctuations with power-spectrum 
P(k) oc k'^T^k), where T{k) = [1 + 6.8k + 72.0k^^^ + 16.0k'^]-^ is the CDM transfer function 
(Davis et al. 1985). More details about these simulations can be found in Moscardini et al. 
(1995). 

The presence of the neutrino component in BMDM was taken into account in the initial 
conditions as follows. We used the transfer functions computed in Section 3 for each component 
separately. At a redshift 2; the neutrinos have a thermal motion with velocity distribution given 
by the Fermi-Dirac statistics 

dn\v) oc 7—. — , (13) 

^ ' exp[v/^;o(;2)] + 1 ^ ' 
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where vo{z) — 7.2(1 + z){'mu/7eV)~^ km s~^, rrii, being the neutrino mass; the corresponding 
rms velocity is t'rmsC-^) — 3.596 Vo{z). We start our simulations ai, z — 20. 

Very rapidly the growth of baryonic fluctuations becomes similar to that of the cold compo- 
nent: because of this reason (see also Klypin et al. 1993) we preferred to include the baryons 
with the cold particles and considered only two components. Our N-body simulation of BMDM 
includes three sets each having 128^ particles, one for the cold particles and two for the hot 
ones; the relative masses of cold and hot particles arc 0.7 and 0.15 respectively. The initial 
conditions were generated by following the same method used by Klypin et al. (1993), still 
based on the Zel'dovich approximation. The displacement of the cold particles was obtained 
with a power-spectrum where only the transfer function of the cold plus baryonic components 
was considered. The same random phases were used to generate the displacement of the hot 
particles: in this case the power-spectrum contains only the transfer function for the hot com- 
ponent. A thermal velocity with random direction and amplitude drawn from the Fermi-Dirac 
distribution was added to the velocity of the hot particles. Imposing that the two sets of hot 
particles have velocities with opposite directions avoids the generation of spurious fluctuations. 

4.1 Matter distribution 

In Figure 4 we plot the projected particle positions from a slice of depth 20 Mpc, for 
the three models at the present time. A 1/4 random sampling is adopted and only the cold 
particles are shown in the case of BMDM. At first glance, the resulting distributions for SCDM 
and BMDM are quite similar, even if in the latter case larger underdense regions are present. 
Structures are instead more evident for BCDM, where voids are larger, the mass distribution 
being more lumpy with many isolated clusters. 

An interesting insight on the clustering of the mass can be obtained from the counts-in-cells 
analysis. In particular, we calculate the variance cr^ = (5^). Figure 5 shows as a function of 
the side R of cubic cells for the three models at the present time; shot-noise corrections were 
applied. For comparison, we also show the results (with 95% error bars) obtained by Loveday 
et al. (1992) in the analysis of the Stromlo-APM galaxy redshift survey. According to Loveday 
et al. (1992) galaxies in this catalog are nearly unbiased with respect to the mass; this implies 
that their counts can be directly related to the matter distribution. The results we show refer 
to real space: we checked that redshift distortions in our simulations slightly affect the variance 
only on very small scales. 

The SCDM model, as already known (Efstathiou et al. 1990; Loveday et al. 1992), has too 
small a variance on scales ranging from 20 to 60 Mpc. The clustering for BCDM is too 
strong on small scales: is outside the 95% observational range for 35 Mpc. On the 
contrary, the results for BMDM are very close to the observational ones. 

One of the most important features of the large-scale matter distribution is the presence 
of large voids: for example, de Lapparent, Geller, & Huchra (1986) found in slices of the CfA 
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survey empty regions of radius around 50 h"^ Mpc. More recently, Vogeley et al. (1994a) 
compared void statistics of the extended CfA survey with simulations of different cosmological 
models, finding that models with Gaussian initial fluctuations and CDM-like primordial power- 
spectra are favoured. Comparing the Perseus-Pisces survey with N-body simulations, Ghigna 
et al. (1994) find that a MDM model with Qcdm/^u/^b = 0.6/0.3/0.1 and primordial spectral 
index n = 1, exceeds the observational Void Probability Function for radii smaller than 10 h^^ 
Mpc, while a CDM model fares better. Statistics aiming at measuring the presence of voids 
in the universe are plagued by the very ambiguity in the definition of 'void'. For our analysis 
we preferred to calculate the Underdensity Probabihty Function (UPF) Pso{R), defined as the 
probability of having a cubic region of side R with a density more than 80% below the mean 
(i.e. S < —0.8). This statistic is less sensitive to the estimate of the mean density and to shot- 
noise, and permits to probe larger scales (being non-zero for larger R) than the usual Void 
Probability Function. Figure 6 shows Pgo as a function of R. As expected, and in agreement 
with the analysis of Piran et al. (1993), a blue primordial spectrum increases the presence of 
large voids. In fact, while the probability for SCDM goes to zero for i? 25 Mpc, BMDM 
and BCDM have a non-zero Pgo out to 30 and 40 h"^ Mpc respectively. For small-scale voids, 
SCDM and BMDM behave in a very similar way; on the contrary, as also evident from the slice 
in Figure 4, BCDM presents also a large amount of small voids. 

Another test of the matter distribution we consider is the mean contour genus per unit 
volume gs- Recently, Vogeley et al. (1994b), performing a topological analysis of the CfA 
redshift survey, found that SCDM has too a large amplitude on smoothing scales ^ 10 
Mpc, while an open CDM model and a CDM model with non-zero cosmological constant are 
always consistent with the observed genus. Figure 7 shows, for the three models at the present 
time, the gs curves obtained using the program by Weinberg (1988). The simulated data 
have been filtered by a Gaussian window, W{r) = (tt^/^A^)"^ exp(— r^/Ag), with two different 
radii, Ae = 6 h"^ and Ae = 20 h'^ Mpc, in order to examine the topology of structure on 
both nonhnear and linear scales. For the ease of comparison with Vogeley et al. (1994b), 
we consider contours of fixed volume fraction and we plot gs as a function of the threshold 
ly parameterizing the volume fraction in a Gaussian field: F^oi = (27r)~^/^ exp(— t^/2)o?i. 
Moreover we normalize the genus curves using the volume of the CfAl+2 survey. Error bars, 
obtained by a bootstrap technique, are always smaller than 10%; for clearness they are not 
shown in the figure. For a Gaussian random field (e.g. Doroshkevich 1970), 

55 = 5o(l-i^')exp(-i/V2), (14) 

where the normalization gg is related to the power-spectrum of density fiuctuations P{k) 
through 

_ 1 
^°-4^ 
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The action of gravity changes the underlying statistics for the density fluctuation 5, breaking 
the symmetry oi gs- A shift of the peak of the genus curve towards either high or low density 
reflects a topology that is bubblc-likc or mcatball-like respectively. In order to quantify this 
shift, we use a statistics Au (e.g. Park, Gott, & da Costa 1992; Vogeley et al. 1994b) defined 
as 

Au = ^V""^^"^"- , (16) 

where g{i')obs and g{iy') fu are the measured and the best-fit random-phase genus curves respec- 
tively. In agreement with the previous definition, Av is positive (negative) when the genus 
curve has a bubble (meatball) shift. In Table 4 we compare the results for go and Au at the 
two filtering radii for our three models and the CfAl-|-2 sample. In agreement with Vogeley et 
al. (1994b), SCDM is found not to reproduce the observed topology in the nonlinear regime: 
the amplitude of go at Ag = 6 h^^ Mpc is too large and the bubble shift too small; moreover, at 
Ae = 20 h^^ Mpc the peak shift is in the wrong direction. A similar negative trend is shown by 
BCDM, both at small and large smoothing. On the contrary, BMDM has the desired behaviour: 
the CfAl-|-2 data are well inside the error bars. 



Table 4. Topological statistics. 





Ae = 6 h-^ Mpc 


Ae = 20h Mpc 




9o 


A, 


9o 


A. 


CfAl+2 
SCDM 
BCDM 
BMDM 


7.70 
11.37 ±2.09 
11.11 ±1.97 
8.47 ± 1.75 


0.19 
0.07 ±0.12 
0.05 ±0.13 
0.13 ±0.12 


4.00 
3.98 ± 1.30 
5.02 ± 1.28 
3.77 ± 1.09 


0.14 
-0.13 ±0.14 
-0.07 ±0.14 
0.05 ±0.16 



4.2 Velocity field statistics 

In order to recover the velocity field from each simulation, we follow a standard procedure 
(e.g. Kofman et al. 1994). First, we interpolate the mass and momentum from the particle 
distribution onto a cubic grid with 128^ grid-points, using a Triangular Shaped Cloud algorithm 
(see e.g. Hockney & Eastwood 1981); we then further smooth by a Gaussian filter with width 
2 h-^ Mpc, to ensure a non-zero density at every grid-point. The velocity at each grid-point 
is defined as the momentum divided by the mass density. 

The left panel of Figure 8 shows the velocity modulus v = \v\ distribution function P{v) 
for our models. The blue models generally produce larger flows. The difference is statistically 
significant: in BCDM and BMDM there are respectively about 700 and 13,000 grid-points with 
a velocity greater than 2000 km s~^, whereas SCDM contains no grid-point with such a high 
velocity. 
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To further characterize the velocity field of each model, we computed the bulk fiow, i.e. 
the amplitude of the center-of-mass velocity of a sphere: Vbuik{R) — \J2i=i'i^i \ /n, where the 
sum extends over the n grid-points falling within a distance R from the center. In order to 
compare our results with observational data from the POTENT analysis (e.g. Dekel 1994) we 
smooth the simulated velocity field by a Gaussian filter of radius 12 Mpc. In particular, 
the bulk fiow is calculated for each model by randomly selecting 100 different grid-points from 
the simulations and calculating the statistics in top-hat spheres, of radius ranging from 10 to 
130 Mpc, centered on them. The central panel of Figure 8 shows the values obtained by 
averaging the results over the 100 estimates. Error bars, which are in all cases less than 5%, 
have not been plotted for clearness. Due to the large error bars on POTENT data, all three 
models are in good agreement with observations. As expected, the larger flows in blue models 
displayed by the P{v) analysis imply higher bulk flows, in better agreement with the POTENT 
central values. This is the reason of the higher likelihood resulting for blue models in the linear 
theory analysis reported in Section 3.1. Of course, having larger data samples, such as the 
"Mark 111" compilation (Willick et al. 1994), can help to increase the discriminatory power of 
this statistical test on the large-scale velocity held. 

Another relevant statistic for the peculiar velocity fleld is the Cosmic Mach Number, Ai 
(Ostriker & Suto 1990), deflned as the ratio of the bulk flow to the one-point velocity dispersion 
(7^: Ai{FC} — VbuikiR) / CTviR) , where cr^(i?) = ^i=iivi — viy^ikf /n and the sum is extended over 
the same grid-points. Since the bulk flow is caused by density fluctuations on scales larger 
than the sampled volume, while the velocity dispersion mostly depends on the power on smaller 
scales, M. actually measures the ratio of large to small-scale power in the velocity fleld. The 
mean values of M., obtained in the same way as for the bulk flows, are also shown in Figure 
8 (right panel), for the three models at the present time (where error bars, always less than 
10%, are not plotted). The different estimates are very close and always inside the error bars. 
This result shows that the Cosmic Mach Number poorly discriminates among different models: 
a similar conclusion was obtained in the study of non-scale-invariant (Moscardini et al. 1995) 
and non-Gaussian (Lucchin et al. 1995) GDM models. 

4.3 Cluster analysis 

In order to identify clusters in our simulations we adopted a standard method (White et al. 
1987; Jing et al. 1994; Jing & Fang 1994; Groft & Efstathiou 1994). As a first step, we adopt 
a friends-of-friends algorithm to extract the groups of cold particles with a linking parameter 
equal to one fifth of the mean separation of the cold particles. After finding the center of 
mass of the groups, we calculate the total (i.e. cold plus hot) mass inside a sphere of radius 
Tci = 1.5 h^^ Mpc, roughly corresponding to the radius of Abell clusters. When the distance 
between the centers of two cluster candidates is smaller than 2rc/, we eliminate from the list 
the cluster with the smaller mass and we repeat the whole procedure using all particles not 
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belonging to 'clusters', until the number of clusters does not change: in this way we can take 
into account the possible presence of more clusters in overdense regions and/or filaments (i.e. 
superclusters) which would be merged by the percolation algorithm. 

In Figure 9 we show the spatial distribution of clusters obtained in our simulations. In 
this plot, we consider only the richest clusters, i.e. those with a mass larger than a suitable 
threshold chosen to reproduce the mean density of the ACO/Abell catalogue, corresponding to 
an intercluster separation dc ~ 40 Mpc. 

The cumulative mass function n(> M), i.e. the number of clusters per unit volume with a 
mass larger than M, is presented for the three models in Figure 10. The BCDM model has a 
larger number of rich clusters compared with SCDM and BMDM, which behave in a similar way. 
A comparison with the observational estimates can be done using the mass fTinction derived 
by Bahcall & Cen (1992) using both optical (richness, velocities and the luminosity function of 
galaxies in clusters) and X-ray data (temperature distribution of clusters). All models present 
too a high mass function for larger masses. The disagreement, already known in the case of 
SCDM (e.g. Bahcall & Cen 1992; Jing & Fang 1994) is significant: for example, no object 
with mass larger than 10^^ Mq is found in the Bahcall & Cen (1992) data, while these objects 
may be found in the simulations. However, the estimates of Bahcall & Cen (1992) should be 
reconsidered, because of some underlying assumptions, such as the density profile or sphericity 
(e.g. Klypin & Rhee 1994). Independent estimates of the cluster mass function, using masses 
derived directly from the dynamics (i.e. virial evaluations), give lower values for the mass 
function (Biviano et al. 1993). Data from Bahcall & Cen (1992) and Biviano et al. (1993) 
are also shown in Figure 10; in the latter case error bars account for uncertainties both in the 
best-fit parameters and in the cluster mean density. Considering the Biviano et al. (1993) 
data, the agreement with BMDM and SCDM is definitely better, while the number of very rich 
cluster is still a problem for BCDM model. 

These results are in good agreement with the linear cluster analysis discussed in Section 
3.2. All the considered models tend to overestimate the cluster abundance, as it is clear from 
the fact that the Press-Schechter formula fits the Bahcall & Cen (1992) data only provided the 
threshold is chosen to take high values, such as Sc ~ 2.5, which is equivalent to an artificially 
low normalization (high b) for these models. If, on the other hand, the cluster abundance had 
been used as a constraint to fix the power-spectrum normalization (as done by e.g. Pogosyan & 
Starobinsky 1995), a value of erg in contrast with the amplitude required by COBE data would 
have come out. 

For our simulated cluster catalogs we compute also the two-point correlation function ^cc- 
Observational estimates of ^cc are well fitted by the relation — {f/fo}"^'^- The amphtude 
of the correlation function and its possible dependence on the cluster richness are rather con- 
troversial: recent estimates of tq based on different samples range from 13 to 25 h'^ Mpc (e.g. 
Nichol, Briel, & Henry 1994 and references therein). Figure 11 shows the results for the cluster 
catalogs in our simulations. For clearness, error bars, obtained by bootstrap resampling, are 
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only shown for BMDM. As already known (see e.g. Bahcall & Cen 1992), the cluster correlation 
function for SCDM is too low. On the other hand, BCDM cannot reproduce the right slope: in 
fact, this model is overclustered at small scales and underclustered at larger scales. Only the 
results for BMDM are well inside the observational range: in this case ro is approximately 16 
h-^ Mpc. 

5 Conclusions 

In this work we presented a detailed comparison of cold and mixed dark matter models, with 
blue primordial power-spectra, with observational data on large scales. In our analysis we used 
both linear and nonlinear techniques. Let us stress that this new class of models, based on 
the assumption that n > 1, has the same level of naturalness as those based on any other 
choice of the primordial spectral index, such as the usual scale-invariant n — 1 value or the 
tilted, n < 1, ones. Being primordially scalc-frcc our initial conditions are more natural than 
those where the power-spectrum contains a built-in characteristic scale. Blue types of spectra 
are suggested by the need of reconciling the amount of power required on large (i.e. in the 
range 10 — 100 Mpc) scales with the amplitude of perturbations implied by COBE data on 
larger scales. An unavoidable consequence of this choice, in the frame of cold dark matter 
models, is an excessive level of fluctuations on small (i.e. ^ 1 Mpc) scale. This is the strongest 
motivation for resorting to the free-streaming effect of a hot dark matter component. Indeed, 
the models which fare better in most of our analyses are the blue mixed dark matter (BMDM) 
ones. Since our analysis assumes the COBE DMR data as a normalization constraint, the first 
test to consider is the level of temperature fluctuations on degree scale. De Bernardis, Balbi, 

6 Vittorio (1995) tested BMDM models with a late, sudden reionization of the intergalactic 
medium against the ARGO (de Bernardis et al. 1994), COBE DMR (Bennett et al. 1994), 
MAX (Devhn et al. 1994) and MSAM (Cheng et al. 1994) data set. Their likelihood analysis 
indicates that mixed dark matter models with blue power spectra (n ~ 1.2) and a reionization 
at z ~ 30 are most consistent with the data, even if the standard model with n — 1 and without 
late reionization is statistically indistinguishable from the best model. It is not a problem to 
have a late reionization in the context of blue models. Indeed, they can have enough small-scale 
power at early times to produce an early generation of structures, leading to late reionization 
(e.g. Carr, Gilbert, & Lidsey 1994). Related to this property is also an anticipation of the 
galaxy and cluster formation epochs compared to the scale-invariant MDM case (see, e.g., 
Lyth & Liddle 1995). 

The large-scale matter distribution emerging from our BMDM scenario has the right prop- 
erties to fit the existing observational data on various statistics, such as the counts-in-cells 
analysis of Stromlo-APM galaxies and the topology of the CfA redshift survey; thanks to the 
increased power on large scales compared to SCDM and MDM, the large-scale voids in our 
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simulations also appear in good shape to reproduce the amount and size of those observed in 
galaxy catalogs. We also computed, both within linear theory and by numerical simulations, 

the large-scale peculiar velocity field, and performed a statistical comparison with observa- 
tional data from POTENT (e.g. Dekel 1994): although the best model selected by a Likelihood 
analysis is the blue CDM one, wc find that BMDM models, with a small fraction of hot dark 
matter can generally provide a reasonable agreement with the peculiar velocity data. 

A problem which remains unsolved in the frame of all the models we considered is related 
to the cluster abundance. The general tendency is to produce too many massive clusters. 
Our analysis, in this respect, strongly differs from that of Pogosyan & Starobinsky (1995), 
who normalized their power-spectra to fit the cluster abundance, as deduced by a Press- 
Schcchter formula, with the Bahcall & Cen (1992) data. For us, instead, once the COBE 
DMR normalization has been imposed on the models, we obtain the cluster mass function as 
a prediction of our analysis. We believe, on the other hand, that many aspects of the cluster 
distribution represent an open problem for all natural (i.e. not containing ad hoc parameters) 
scenarios of structure formation, whose solution would probably require a better modeling of 
cluster formation in numerical simulations than it is presently possible by coUisionless codes. 
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Figure captions 



Figure 1. The total transfer function (sliort-daslied line) and the transfer function for bary- 
onic (continuous line), cold (dotted line) and hot components (long-dashed line), at different 
redshifts: a) z — 20; h) z = 10; c) the present time z — 0. The results refer to a fiat MDM 
model with — 0.3, Qb — 0.05 and h — 0.5. 

Figure 2. The rms value of the fluctuation density field at the present epoch, a{R), is plotted 
versus the filtering scale R, in units of Mpc {h — 0.5 has been used). The four curves refer 
to the different scenarios of galaxy formation that we consider in this paper: BCDM (dotted 
curve), SCDM (dot-dashed curve), BMDM (continuous curve) and MDM (dashed curve). All 
these models have initial spectra normalized to the rms temperature fluctuation of the COBE 
DMR maps. 

Figure 3. The predictions for the local X-ray luminosity function are compared with data 
from Kowalski et al. (1984). The dotted curves in the three panels arc evaluated with the 
parameter choice listed in Table 3 which give the best-fit of the data (sec text for details). The 
dashed curves in each panel arc the XRLFs obtained with 5c = 1.7 and fluctuation spectra 
normalized to the COBE/DMR amplitude. X-ray luminosities are evaluated in the 2 — 6 keV 
energy band. 

Figure 4. Matter distribution in shces of 20 h'^ Mpc for SCDM (top panel), BMDM (left 
down) and BCDM (right down), at the present time. The box refer to the whole simulation, 
260 Mpc. A random sampling of 1/4 is adopted. For BMDM only the "cold" particles are 
shown. 

Figure 5. Counts-in-cells analysis: the variance as a function of the side R of cubic cells 
for SCDM (sohd line), BMDM (dotted line) and BCDM (dashed hne). The filled circles and 
the heavy sohd lines refer to the Loveday et al. (1992) analysis of the APM-Stromlo catalog 
and their 95% confidence level, respectively. 

Figure 6. The Underdensity Probabhty Function Pgo as a function of the side R of cubic cells 
for SCDM (solid hne), BMDM (dotted line) and BCDM (dashed line). 

Figure 7. The mean contour genus as a function of the effective threshold u for SCDM (solid 
hne), BMDM (dotted hne) and BCDM (dashed hne). Two different Gaussian filtering radii are 
considered: Ae = 6 Mpc (left panel) and Ae = 20 Mpc (right panel). 

Figure 8. Velocity field statistics. Different lines refer to different models: SCDM (sohd line), 
BMDM (dotted hne) and BCDM (dashed line). Left panel: the distribution function of the 
velocity modulus v. Central panel: the bulk flow as a function of the radius R of the sphere. 
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Filled circles and error bars refer to the POTENT data (Dekel 1994). Right panel: the Cosmic 
Mach Number as a function of the radius it! of the sphere. 

Figure 9. Projected distribution of clusters identified with an intercluster distance of 40 
Mpc. The box represent the whole simulation (260 h'^ Mpc side): SCDM (left panel), BMDM 
(central panel) and BCDM (right panel). 

Figure 10. The cluster mass function N > (M) (M in units of the solar mass) for SCDM 
(solid line), BMDM (dotted line) and BCDM (dashed line). For comparison the observational 
estimates by Bahcall & Cen (1992) and Biviano et al. (1993) are shown by open circles and 
filled squares, respectively. 

Figure 11. The cluster correlation function ^cc as a function of the distance r for SCDM (solid 
line and open squares), BMDM (dotted line and open circles) and BCDM (dashed line and 
filled circles). Bootstrap errors are shown only for BMDM. As comparison, the two heavy solid 
lines represent the power-law ^cc = {1^/1^0)'^'^, with tq = 13 and 25 Mpc. 
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